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SUBLIMATION  OP  A  BLUNT  BODY  IN  THE  VICINITY  OF  THE 
STAGNATION  POINT  IN  PLANE  AND  AXI SYMMETRICAL  FLOW 
OF  A  GAS  MIXTURE 

G.  A.  Tirskiy 
(Moscow) 

The  problem  of  the  sublimation  of  a  body  In  the  vicinity  of  the 
stagnation  point  was  examined  in  an  earlier  work  [l],  where  a  numerical 
solution  was  carried  out  for  the  case  of  quasistationary  "adiabatic" 
(disregarding  heat  transfer  into  the  body)  boiling  of  solid  carbon 
dioxide  in  an axi symmetrical  case  where  the  Mach  number  of  the  incident 
flow  was  6.2.  However,  the  conditions  for  attaining  such  a  maximum 
regime  of  sublimation  were  not  obtained  in  the  work  and  no  solution 
was  given  which  permitted  the  calculation  of  the  sublimation  of  bodies 
with  other  physical  properties  under  arbitrary  flow  conditions. 

A  simplified  analysis  of  the  sublimation  of  a  solid  at  a  given 
sublimation  temperature  T^  has  been  presented  [2]  which  is  based  on 
the  method  of  integral  relationships.  This  case  also  corresponds  to 
the  "boiling"  of  a  solid  since  during  sublimation,  the  temperature 
prior  to  boiling  temperature  on  the  vaporizing  surface  is  unknown  and 
must  be  determined  while  solving  the  problem. 
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In  the  present  article  we  obtained  an  accurate  solution  (solution 
of  the  type  of  a  uniformly  propagating  wave  )  of  the  problem  of 
equilibrium  and  nonequilibrium  sublimation  of  a  blunt  body  in  the 
stagnation  region  with  an  arbitrary  dependence  of  the  physical  prop¬ 
erties  of  the  body  on  tenperature. 

For  a  Prandtl  number  a  ■  0.7;  1  and  a  Schmidt  number  Sc  -  1  the 
solution  was  obtained  numerically.  By  means  of  asynptotic  integra¬ 
tion  a  solution  was  also  obtained  for  arbitrary  a  and  Sc  >  0.5.  It 
was  demonstrated  that  if  the  coefficient  of  accomodation  f  >  0.1  then, 
with  sufficient  accuracy  for  application,  vaporization  of  the  body  will 
proceed  according  to  diffusion  kinetics  (equilibrium  vaporization) 
for  ~  10s  m/sec,  when  f  <  0.1  nonequilibrium  vaporization  must  be 
taken  into  account. 

The  necessary  and  sufficient  condition  under  which  boiling  at  the 
sublimation  front  takes  place  was  obtained  (7) .  When  boiling  Is 
attained  the  sublimation  rate  and  the  mass  velocity  of  sublimation  are 
In  a  final  form  (Eqs.  7.4  and  7.5). 

The  temperature  profile  in  the  body  for  an  arbitrary  dependence 
of  thermophysical  properties  of  the  body  on  temperature  are  found  in 
quadratures.  The  calculation  of  any  specific  problem  in  the  general 
case  reduces  to  a  solution  of  a  system  of  three  finite  equations  for 
the  determination  of  concentration,  temperature  at  the  vaporization 
front,  and  rate  of  sublimation. 

1.  Statement  of  the  Problem 

If  a  mixture  of  gas  flows  past  a  solid  body  with  a  vapor  partial 
pressure  of  the  body  in  the  Incident  flow  lower  than  the  pressure  of 
the  saturated  vapors  at  the  body  surface  temperature  and  if  the  sur¬ 
face  temperature  is  lower  than  the  temperature  of  the  triple -point 
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of  the  phase  diagram,  then  the  body  will  sublimate.  For  this  purpose 
we  will  work  out  for  the  sublimation  front,  the  appropriate  tempera¬ 
ture  and  mixture  composition  which,  like  the  position  of  the  sublima¬ 
tion  front,  are  not  known  beforehand  and  must  be  determined  while 
solving  the  problem. 

The  process  of  sublimation  of  a  multicomponent  body  must  be  con¬ 
sidered  as  a  simple  heterogeneous  chemical  reaction  with  the  stoichio¬ 
metric  equation 


V.  v'k  a\ 

km  1 


►  2  v* 
tm\ 


(1.1) 


where  A^,  are  the  chemical  symbols  corresponding  to  the  solid  and 
gaseous  components;  v^,  v^  are  the  stoichiometric  coefficients  of  the 
solid  and  gaseous  components  respectively;  N0  13  the  number  of  com¬ 
ponents  In  the  solid  phase;  and  N‘  Is  the  number  of  .  v-.ponent3  pro¬ 
duced  during  sublimation  of  the  body  (for  example,  during  the  vaporiza¬ 
tion  of  S102  a  four -component  ion  vapor  of  Si02,  SiO+,  ot  and  0+  is 
formed  [3]) . 

The  mass  concentrations  c^s(i  =  1>  . ..,  N«)  of  the  gaseous  pro¬ 
ducts  approaching  the  sublimation  surface  from  the  side  of  the  solid 
phase  are  associated  with  N*  by  the  equation 


v».V, 


cAr', 


\N.  Mn.  ’  — ^ 


V,  Ck,  =  1, 


where  is  the  molar  mass  of  the  i-th  component.  Hence  the  composi¬ 
tion  of  the  gaseous  products  produced  by  vaporization  of  the  body  will 


be 


Ci.  = 


I  vt.u4.' 


i . y. 


(1.2) 


Then  the  problem  of  a  stationary  regime  of  sublimation  in  the  stagna¬ 
tion  line  (plane  case)  or  in  the  stagnation  region  (axisymmetrical  case) 
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of  a  blunt  body  placed  In  a  stationary  gas  flow  consisting  of  N  »  N*  + 


+  N"  components  leads  to  the  solution  of  a  system  of  equations  of  a 
nonstationary  boundary  layer  for  a  multicomponent  gas  mixture: 


!•;  -r.^T  ^  (P"*n_1)  t  Ty  (?»’>  =  0. 
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(1.5) 

(1.6) 
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together  with  the  equation  of  thermal  conductivity  in  the  solid  body 


f V.?T -£(*«?£)’  (1.8) 

Here  x,  jr  are  the  coordinates  associated  with  the  body  along  and 

normal  to  the  initial  surface  of  the  body  respectively;  u  and  v  are 

the  components  of  the  velocity  vector  with  respect  to  x  and  jr 

respectively;  p  ia  the  density;  p  -  pressure;  T  and  Ti  —  temperature  of 

the  gas  and. body;  c^  —  mass  concentration  of  the  i-th  component  of 

the  mixture;  h,  .=  h?  +  h?  =  c_  T  +  h?;  and  c  —  specific  partial 
.1  1  l  p^  i  p^ 

enthropy  and  specific  heat  for  constant  pressure  of  the  i-th  component 
respectively;  h°  —  "zero  enthalpy"  of  the  i-th  component  [4];  p.,  X, 

Diz  —  coefficients  of  viscosity,  thermal  conductivity  and  diffusion 
of  the  mixture;  R  —  universal  gas  constant;  pi,  c#,  Xi  —  density, 
heat  capacity,  and  coefficient  of  thermal  conductivity  of  the  solid 
body;  L(Ti)  and  N(Ti)  —  given  functions  of  temperature. 
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It  is  easy  to  show  that  in  an  approximation  of  the  boundary-layer 
theory  a  term  with  baro  diffusion  drops  out  of  the  equations  of  dif¬ 
fusion  (1.5)  and  the  equation  of  energy  (1.6). 

When  writing  equations  (1.5) -(1.7)  we  must  also  assume,  in  addi- 
'  tion  to  the  usual  assumptions  of  boundary-layer  theory  in  the  stagna¬ 
tion  region  that  all  of  the  Maxwellian  diffusion  coefficients  are 
equal  (for  a  binary  mixture  this  supposition  is  satisfied  exactly) , 
that  the  effect  of  thermal  diffusion  is  small,  and  that  the  mixture  is 
ideal.  The  equation  of  thermal  conductivity  (1.8)  is  written  on  the 
assumption  that  the  thickness  of  the  thermal  boundary  layer  in  the 
body  is  much  smaller  than  the  radius  of  curvature  of  the  surface  of 
the  body  in  the  vicinity  of  the  stagnation  point  (line) . 

System  of  Equations  (l.3)-(l.8)  is  solved  for  the  following 
boundary-value  conditions: 

1)  on  the  outer  edge  of  the  boundary  layer 

U  =  ?X,  Ct  =  Ci»  (i«l,  ...  ..V),  h=  ha,;’-  (1.9) 

2.  from  the  law  of  conservation  of  mass  find  energy  on  the  surface 
of  sublimation  which  is  unknown  before  the  solution  of  the  problem 
we  obtain  [15]  p  (D  —  v)  *  pxD; 

dr  1  _ 

pci  (D  —  p)  -f  pDn  =  PiCuD  (i  = cl... ..An.  (l.io) 

•  ptj  {O  r)  pDia  ~  =  0  (/ <*  1. ... ,  .V); 

•  »<»-•»[«  <« + -s  <-*•)]  p d„  l.  a 

or 

p  (D  ~  v)  [l  (T9)  +  Lc  (*•  -  V  /*:)]  =  PD»  Le  ^ ; 

*”* 

3)  at  infinity,  in  the  solid  body 

r,  -r_,  (1.11) 
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where  D  is  the  displacement  velocity ,  normal  to  the  surface,  of  the 
sublimation  front  relative  to  the  body  and  l (To)  is  the  latent  heat 
of  sublimation  which  is  dependent,  generally  speaking,  on  the  surface 
of  sublimation  To. 

Since  hereafter  we  will  trace  the  stationary  regime  of  sublima¬ 
tion,  we  will  not  put  forth  the  initial  conditions.  The  formulated 
problem  is  a  parabolic  system  of  nonlinear  equations  (1.3) -(1.6) , 

(1.8)  of  the  2N  +  7th  order  with  2N  +  5  conditions  (l.9)-(l.H)  which 
contain,  in  addition,  the  two  unknown  quantities  T0  and  D.  In  order 
to  complete  the  problem  it  is  necessary  to  introduce  four  additional 
relationships  which  characterize  the  kinetics  of  sublimation. 

For  a  sufficiently  dense  viscous  mixture  we  can,  of  course, 
assume  an  equality  of  tangential  velocities  at  the  sublimation  fronti 
i.e.  . 

U  -  0.  (1.12) 

It  is  known  [6]  that  during  weak  vaporization,  I.e.,  when  con¬ 
vection  heat  transfer  from  the  vaporizing  surface  Is  significantly 

•  '  /•  .  1  * 

lower  than  molecular  heat  transfer,  the  temperature  Jump  Is  of  the 
order  of  X/a  (the  ratio  of  the  length  of  free  path  to  the  character¬ 
istic  dimension  of  the  problem),  i.e.,  T  =  Ti  =  TQ.  When  convection 
heat  transfer  is.  significantly  greater  than  molecular  heat  transfer, 
the  temperature  Jump  and  deviation  of  vapor  partial  pressure.  p0  form 
equilibrium  p(o)'  at  the  sublimation  front  will  be  of  the  order  of  v/c 
(the  ratio  of  the  vapor  flow  velocity  to  the  average  velocity  of  the 
thermal  motion  of  the  molecules) .  If  it  is  assumed  that  during 
vaporization  the  thickness  of  the  boundary  layer  remains  small  with 
respect  to  the  characteristic  linear  dimension,  then  v  ~~  (?  —  ujd), 

where  u^  is  the  velocity  of  an  undisturbed  flow,  d  is  the  radius  of 
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curvature  of  blunting,  and  v  is  the  kinematic  coefficient  of  viscosity. 
Then  v/c  «  1,  (c  ~  10*  to  10*  m/sec)  and  in  this  case  we  can  assume 

T  -  Ti  -  T0  (1.13) 


and 

/>„  =  />«»  =/>,-•••+  Ps-  -  W[r)' 


where  the  last  relationship  is  the  equation  of  the  vapor  tension  curve; 
p#  and  T#  are  some  values  of  pressure  and  temper ature  on  this  curve; 
p^(i  “1,  N*)  is  the  partial  pressure  of  the  i-th  vapor  component. 

If  Jhass  concentrations  are  introduced  then  the  equation  of  the  vapor 
tension  curve  may  be  written  in  the  form 


(1.1*) 


For  a  binary  system  .1,  J.,  with  use  of  the  van  der  Vaals  formula  [4] 
for  the  vapor  tension  curve,  Eq.  (1.14)  is  rewritten  as 


(1.15) 


In  order  to  clarify  the  Influence  of  nonequilibrium  vaporization,  I.e., 
when  p0  ^  p(o),  in  place  of  (1.15)  we  will  use  the  fundamental  formula 
of  Knudsen-Langmulr  for  the  mass  velocity  of  vaporization: 


=  /  j/" l/>„  -  p  (0)).  p(0)  =  pjp  ( ,  (l  .16) 

where  f  Is  the  coefficient  of  accomodation  expressing  probability  that 
a  molecule  of  vapor  striking  the  vaporization  surface  will  adhere  to  it 
Formula  (l.l6)  was  derived  on  the  assumption  that  the  flow  of  mol¬ 
ecules  striking  vaporizing  surface  has  a  maximum  distribution  for  a 
quiescent  gas  with  a  temperature  T0. 

As  has  been  demonstrated  [6],  for  a  mass  vaporization  rate  (if 


the  temperature  Jump  Is  neglected)  when  using  the  distribution  function 
in  an  approximation  of  "13  moments"  for  the  calculation  of  the  molec¬ 
ular  flux  impenging  on  the  surface,  a  formula  analogous  to  (l.l6)  is 
obtained  but  in  place  of  £,  we  must  use  2f.  Since  the  law  according 
to  which  dispersion  of  the  molecules  of  vapor  takes  place  is  unknown, 
we  will  use  the  relationship  in  (l.l6).  We  note  that  certain  authors 
[7]  maintain  that  f  ~  1..  Introducing  the  molar  vapor  .concentration 
c*  into  (l.l6)  we  obtain  for  a.  binary  gas -vapor  mixture 


V  -  //•  (0)  (/f  -  1)  j/ 


r~w 


0 

co 

e*(0) 


(1.17) 


Henceforth  we  will  consider  nonequilibrium  vaporization  for  the  case 
of  a  binary  gas -vapor  mixture  and  use  Eq.  (I.17).  Thus,  additional 
conditions  (1.12),  (1*13),  (1.14)  or,  for  the  case  of  nonequilibrium 
vaporization,  (1.17)  close  the  problem.  The  effect  of  dissociation 
and  radiation  are  not  taken  into  account  when  formulating  the  problem. 
It  Is  known  [8]  that  for  flow  stagnation  in  the  vicinity  of  the 
stagnation  point,  the  influence,  of  equilibrium  dissociation  and  radia¬ 
tion  on  the  magnitude  of  the  thermal  flux  is  small  If  the  temperature 
of  the  wall  Is  lower  than  the  temperature  of  dissociation  (l600-2200oK) 
But  at  such  wall  temperatures  sublimation  for  most  materials  may  attain 
a  significant  value. 


2.  Coefficient  of  Molecular  Transfer  and  Prandt  and  Lewis  Numbers 

Generally  speaking,  accurate  expressions  for  the  coefficients  of 
transfer  (coefficients  of  viscosity,  thermal  conductivity,  diffusion, 
and  thermal  diffusion)  for  gases  and  gas  mixtures  may  be  obtained  by 
calculations  according  to  the  law  of  statistical  mechanics  and  the 
kinetic  theory  of  gases  when  the  necessary  data  concerning  lnter- 
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molecular  forces  Is  available.  However ,  these  formulas  are  too  com- 

•  # 

plex  for  practical  calculations. 

Recently  approximation  formulas  have  been  derived  from  exact 
formulas  for  the  coefficients  of  viscosity  [9],  and  thermal  con¬ 
ductivity  [9*10]  of  •  multicomponent  gas  mixtures,  which  over  a  wide 
range  of  temperatures  [11]  give  values  differing  little  from  the 
exact  formulas  of  the  kinetic .  theory  of  gases.  Conpar is on  of  these 
formulas  with  experiments  [9*10]  also  gives  satisfactory  agreement. 
These  formulas,  in  essence,  are  corrected  first  approximation  obtained 
from  a  rigorous  kinetic  theory  of  gases  and  are  very  convenient  for 
numerical  calculations.  In  order  to  calculate  the  coefficient  of 
viscosity  of  a  mixture  it  Is  necessary  to  know  only  the  molecular 
weight  and  coefficient  of  viscosity  of  the  Individual  components  of 
the  mixture  at  the  same  temperature  at  which  the  gas  mixture  is  found. 
In  order  to  calculate  the  coefficient  of  thermal  conductivity,  we 
need  to  know  the  molecular  weight,  the  coefficients  of  thermal  con¬ 
ductivity  of  the  pure  components,  and  also  either  the  viscosity  of 
the  individual  components  - or' the  specific  heat  of  the  components  at 
the  same  temperature  at  which  the  gas  mixture  is  found. 

1..  Mixture  viscosity.  The  approximation  formula  for  the  coef¬ 
ficient  of  viscosity  of  a  mixture  obtained  from  the  elementary  kinetic 
theory  of  gases  [12]  and  verified  from  considerations  of  the  strict 
kinetic  theory  of  gases  [19]  has  the  form 


Kb  ,-Ui  .  .  (»•)*  y  (2.2;  rj)  /  (t,)  kT 

Kb  “  S  >  x* =  ’  •  T‘ "  ^  ' 


(2.1) 

(2.2) 
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where  -  coefficient  of  viscosity  of  the  individual  components; 

(r0)1  =  (r°)ii  “  distance  between  molecules  at  which  the  energy  of 
interaction  is  zero;  —  absolute  value  of  the  maximum  energy 

Of  attraction;  k  —  Boltzmann  constant;  Y(2;  2;  r^)  —  average  reduced 
collision  cross-section  which  depends  on  the  temperature  and  potential 
energy  of  molecular  interaction  i(t^)  —  correcting  multiplier,  depend¬ 
ent  on  T£,  which,  for  a  wide  temperature  range,  differs  from  unity  .. 
by  less  than  +  0.08  [13]. 

When  T  •  >  5  the  function  Y  (2;  2;  T;)  is  slightly  dependent  on. 
temperature.  Therefore  the  magnitude  of  .x^,  within  the  limits  of 
accuracy  of  formula  (2.1)  may  be  considered  constant  for  the  given 
vapor .  Actually,  for  example,  for  vapors  of  car-bon  dioxide  C02  (1) 
and  air  (j)  we  have  [13]  [(r0)±  =  3.996a,  (r0j j  =  3.617,  Ej/k  =  190°K 
and  Ey/:  »  97°K]: 

=  0.S55  (500tK),  =  0,892  (1000°K),  =  0,904  (2000°K), 

...  =  0f907  (4000°K).. 

Whence  Xj^  z  0.729  for  a  mixture  of  CQ2  and  air  in  this  temperature 
range  (disregarding  dissociation) . 

Formula  (2.1)  may  then  be  written  as  the  product  of  the  function 
which  depends  only  on  temperature  times  the  function  which  depends 
only  on  compositions 

.u  =  Hi  VO  (a);  (2.3) 


where  p^  (T)  is  the  coefficient  of  viscosity  of  some  component  of  the 
mixture,  -  . 

,,M  _  v  (,  4.  f  ^  A)-.  -g  +  j,  c„*y* .  (2.4) 

*— I  j*i  *  (-1  3  jii  1 
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2.  For  the  coefficient  of  thermal  conductivity  of  the  mixture 
we  will  use  an  approximation  formula  derived  from  the  expressions  of 
the  rigorous  kinetic  theory  of  gases  [10].  For  monoatomic  gases  or 
for  gas  mixtures  with  frozen  Internal  degrees  of  freedom  this  formula 
has  the  form 


■-i  j+i  \ 

1.065  .  ■  ,  /A? \  (Ut V*  1  V 

G,>  2  V2  ^  +  \.v j)  J/Vj' 


(2.5) 


A  comparison  of  this  formula  with  the  exact  formulas  of  the  kinetic 
theory  of  gases  and  with  experiment  yields  satisfactory  agreement  [10], 
For  a  mixture  of  polyatomic  gases  the  formula  takes  the  form 

•V  V 


i’l  ivi  * 

When  «  X°  Formula  (2.6)  changes  Into  (2.5).  Since 

'i J, 
x?  "  -ViT?)  x,if 

then  from  (2.5)  and  (2.6)  it  follows  that 


(2.6) 


Gij  =  i,065Gij. 


(2.7) 


The  ratio  X^Aj  for  polyatomic  molecules  with  consideration  of  the 
internal  degrees  of  freedom  is  given  by  the  expression 


~  ~  [wj  Ei>~  T.'  -0.115  +  0,354 -y  , 

where  is  a  correction  multiplier  recently  refined  by  Hirschfelder 
[14]  for  gases  with  polyatomic  chemically  reactive  molecules  and 
molecules  with  excited-  electron  levels. 

Formula  (2.6)  may  also  be  written  In  the  form 


.it. 


E, 


n 


X  «=  (f)  fx  fa). 


(2.8) 
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where  X^T)  is  the  coefficient  of  thermal  conductivity  of  some  com¬ 
ponent  of  the  mixture 


*&«*(  l 


j+i 


(2.9) 


3.  The  diffusion  coefficients  in  a  multicomponent  mixture  are 
expressed  by  the  diffusion  coefficients  of  a  binary  mixture  of  all 
possible  pairs  of  components  of  the  mixture,  by  the  molar  concentra¬ 
tions,  and  by  molecular  weights  of  the  components  proportional  to  the 
corresponding  binary  diffusion  coefficients  [15]. 

4 .  Schmidt  number  Sc  =  p./pDi2 .  By  virtue  of  the  structure  of 
Formula  (2.3)  it  suffices  to  calculate  the  complex  p-j/pD^;  then  it 
is  possible  to  calculate  p/pD^j  also.  We  have  . 


i~i  *  , 


V 


V. 


(2.10) 


where 


■A  V 

[r.))  y(2.2;rJ.)  ^  — 


I  ■ 

1  (r0)i  ( ^0) j  1 1  Tjj  =  —  ,  '  €jj  =  E(Ej. 

"  ei; 


position  of  the  mixture  [13]  and  its  value  varies  between  0.97  and  1 

for  a  wide  range  of  temperature  and  composition  changes  [13] «  The 

ratio  Y(l.l j  t^/Y  (2.2;  Tj)  is  also  weakly  dependent  on  tenperature 

for  sufficiently  large  values  of  and  Tj  (r  >  3) •  F°r  example,  for 

a  mixture  of  molecular  hydrogen  and  air  this  ratio  ..  equal  to  0.830 

for  T  =  600°K  and  0.819  for  T  =  2400°K  etc.  Therefore,  the  value  of 

x?.  ,  within  the  limits  of  accuracy  of  Formula  (2.3)  may  be  considered 
1 J  >  J 

constant  and  characteristic  for  the  given  pair  ij .  For  example,  the 
Schmidt  number  for  a  binary  mixture  will  be 
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(2.11) 


Sc  “  iYz  [* "  ei(j/j~  'jjl1  +  Iff)  G,t "f 

+  (*  +ci‘rr7;)  ]• 

i.e.,  It  Will  not  depend  on  temperature  buy  only  on  composition.  The 
maximum  variation  of  Sc  with  respect  to  £  Is  equal  to 

Sc  =»  1J  \Alj) 

for  example,  for  vapors  of  hydrogen  and  air  ~  0.9. 

5.  Lewis  number ,  Le  *  X/pCpD^.  By  virtue  of  the  structure  of 
Formula  (2.8)  it  suffices  to  calculate  the  complex  *-j/PcpDij*  then  it 
Is  posislble  to  calculate  Le  also.  We  have 

•*  .  (*■*> 

The  maximum  variation  of  Le  for  a  binary  mixture  is  equal  to 


.0(^-0)  .  ePi  Ej  .  <£(*?  , 

Jt  1)  cpj  £j  /  xo  —  <■*»  V.Vj/  Jtij, 


I 

u 


(2.15) 


where  c; 


m 

Pi 


Cp^  is  the  molar  heat  capacity  of  the  i-th  component. 


Since  the  first  two  factors  in  (2.15)  are  close  to  1,  the  variation 
of  Le  with  a  change  of  composition  depends  mainly  on  the  ratio  of  the 
molecular  weights  of  the  components. 

6.  Prandtl  number,  a  =  n.Cp/\  <*  Sc/Le. 

5*  Formulation  of  the  Boundary  Value  Problem  for  a  System  of 
Ordinary  Differential  Equations 


.  Proceeding  from  the  form  of  the  boundary-value  conditions  (1.9) 
and  (l.ll)  and  talcing  into  account  that  the  solution  yielding  a 
stationary  regime  of  sublimation  must  be  a  type  of  uniformly  propaga¬ 
tion  wave,  we  will  search  for  a  solution  of  Eqs.  (l.5)-(l.8)  with 
conditions  (l. 9) -(1.17)  In  the  forms 
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for  a  gas -vapor  mixture 


*;«*;= -r£\. 

A  *Aj-i-(Aoo-A*)a(rj),  et  «  e*  +  (cia8 -el0)g,  (tj). 


(5.1) 


for  a  solid  body 

Tt  -  f-oeOjfa,).  TJX  (y  -/)/),  (5.2) 

where  K  is  an  arbitrary  constant,  the  subscript  0  refers  to  unknown 
values  of  the  parameters  at  the  sublimation  front  and  the  subscript 
»  refers  to  conditions  outside  the  boundary  layer. 

The  law  of  motion  of  the  sublimation  front  for  a  stationary  regime 
of  vaporization  will  be 


yo  -  Dt,  (5.5) 


where  D  is  not  known  before  solution  of  the  problem  of  the  displace¬ 
ment  velocity  of  the  sublimation  front. 

Substituting  Eqs.  (3.1) -(5.5)  into  Eqs.  (l.3)-(l.8)  and  con¬ 
ditions  (1.9) -(I.17)  we  obtain  for  determination  of  the  four  functions 
9,  z,  =  g(i  *  1,  ...,  N),  0i  and  the  parameters  D*,  0O  =  T  o/^t 
cio(i  =1,  N)  the  following  nonlinear  boundary-value  problem: 


K  (ojwp')  '  -r  «<P9'  =  9*  —  PooP**, 

M~l  —  V  — 

f  -u  •  v  ^ 

K  (pfi  Sc"1#')'  -f  nyg’  =  0, 

K  {^JV  4-.  (Le"1  -  1)--”^;  ¥]}'  +  "92'  =  0,  J 


(5.4) 
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I 


[L  (0,)  o;r  -f  D •  OjJV  (0,)  o;  =  0.  g  -  D'o ,  I 


0(O|)*=  5^(61) A. 


(3.5) 


9'  (°c)  “  1,  f  (oc)  a  z  (00)  a  1, 


•9 '  (0)  =0,  g  (0) = c  (0)  a  0,  r,  |/  Jj&L-,  n<p  (0)  ■  D\  0t  (0)  -  01#, 

»<P(0)[2(ru)  +Le#(;4-2  v*:)].+of,jr(*:-ik:)*'w  -o. 


or 

*-»  *-1  * 

— (• ,«*4 


(3.6) 


(3.7) 


(3.8) 


6,  (—OO)  a  1. 


(3.9) 


For  nonequilibrium  vaporization  in  the  case  of  a  binary  gas-vapor 
mixture  with  the  use  of  (1.17)  in  place  of  (5.8)  we  obtain 


jr 

'1 


1/  W4"-  __,m  IpW  ,/ 

V  wZ*’"}m-tZVKV 


(3.io) 


here  the  following  designations  are  introduced: 

X  (T’o)  =  /  (^o)  -r  pjp-'F  7^  [<?  (0,»)  -<?(!)], 


(3.11) 


?*  “j,  P1^.’  r*  =  fa  •  P»  ■*  Pl  (l*o)' 

The  system  (3.^) -(3.5)  la  a  system  of  nonlinear  ordinary  differential 
equations  of  the  9th  order  with.  N  +  2  unknown  parameters  c10  (i  »  1, 

N),  D  and  T0#  which  enter  into  conditions  (5.7)#  (3.8)  or  (j.10) . 
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Relationships  (3.6) -(5. 10)  give  exactly  N  +  11  conditions.  Consequently 
it  may  be  expected  that  the  solution  of  problem  (3.4) -(3. 10)  Is 
determined  uniquely. 

After  solution  of  the  problem,  the  mass  velocity  of  sublimation 
is  determined  by  the  formula 


Pl0D  ==  p10  D *  (0)  p0 


1  /f!£fg-g 

V 


_  p®  „» /r\\  1  /  wA  o. 

--55'  <0)-5^rK  i^?v- 


(3.12) 


The  temperature  profile  in  the  body  Is  found  after  this  in  quadratures 

[15] « 

“  “  \  dir-Wy  (3.13) 


Prom  (3.10)  it  is  possible,  after  solution  of  the  problem,  to 
calculate  the  deviation  of  the  vapor  partial  pressure  p0  over  the 
sublimation  surface  form  the  saturation  pressure  p(o)  at  a  tempera¬ 
ture  T0: 


R  -  jk  -  1  +  "<?  (0)  e 


PnM, 

/p(0)  |/  Af, 


(3.14) 


We  will  evaluate  the  order  of  the  parameter  e.  For  Tw  ~  1000°K, 

'  1  cm8/ sec,  p^  ~  10~4  gr/cm3,  0  ~  103sec-1  (with  uw  ~  103m/sec, 
d  ~  1  m) ,  T0  **•  1000°K  and  p(0)  «*  100  mm  Hg  we  obtain  e  ~  10“5/*  f”1. 

Since  |<p(o)  J  <  1  for  moderate  vaporization,  then  R  *  1  +  10~3f"x. 
Consequently  for  f  >  0.1  it  is  possible  to  obtain  with  a  high  accuracy 
the  equilibrium  vaporization  and  to  use  Eq.  (3.8).  For  f  <  0.01  it 
is  necessary  to  calculate  nonequilibrium  vaporization  and  to  use 
Eq.  (3.10).  In  this  case  the  coefficient  of  accomodation  must  be 
known  exactly. 


This  case  Is  Interesting  in  that  a  number  of  conclusions  may  be 
reached  before  complete  solution  of  problem  (3.4) -(3. 10) .  Since  for 
many  mixtures  the  Lewis  number  is  close  to  unity  (2),  these  conclusions 
will  be  of  practical  Interest. 

Actually,  if  sublimation  proceeds  according  to  diffusion  kinetics 
(condition  3.8)  and  Le  ft  1,  then  z(tj)  s  g(rj)  and  determination  of  the 
condition  and  temperature  at  the  sublimation  surface  leads  to  a  solu¬ 
tion  of  the  final  system  of  N  equations 


unknown  clQ(i  «  1,  ...,  N)  and  To. 

For  a  binary  mixture  (cA  -  c,  c^  «  1  -  c,  c8  -  1)  the  concentra¬ 
tions  of  vapor  and  tenqjerature  at  the  sublimation  surface  are  determined 


from  a  system  of  two  equations 


jo  •  foo  (h,0  !‘in) 


co-c°»  .  _ (4.2) 


It  is  interesting  to  note  that  these  results  do  not  depend  On  the 
character  of  the  dependence  of  transfer  coefficients  on  temperature 
and  composition  and  are  associated  only  with  the  assumption  that  Le  « 

-  i.  _ 

It  follows  from  (4.1)  and  (4.2)  that  "boiling"  (2^J  m  ^  ClQ  »  1) 
on  the  surface  of  the  body  in  the  case  of  equilibrium  vaporization  is 
attained  only  in  the  limiting  case  of  an  infinitely  large  thermal  flux 
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from  the  gas.  Prom  (3.7)  and  (3. 10)  it  follows  that  this  fact  does  not 
depend  on  the  assumption  that  Le  -  1  and  is  associated  only  with  the 
assumption  that  sublimation  proceeds  according  to  diffusion  kinetics. 

If  sublimation  is  nonequilibrium  (f  «  1)  and  the  coefficient  of 
accomodation  is  known,  then  when  "boiling’'  is  attained,  a  solution  to 
the  problem  (3.4) -(3.10)  may  be  obtained  in  final  form. 

Actually  for  "boiling"  p(0)  *  poo  Is  the  stagnation  pressure, 
from  which,  if  the  pressure  at  the  stagnation  point  poo  Is  taken  as 
the  characteristic  pressure  p#,  we  will  obtain  T0  -  T#  *  T0o>  where 
Too  is  the  temperature  on  the  vapor  tension  curve  which  corresponds 
to  the  pressure  of  the  vapors  p0o.  Prom  the  finst  equation  in  (4.2), 
substituting  in  the  left  hand  part  T0  -  T00,  we  find  the  concentration 
of  vapors  cQ  on  the  sublimation  surface.  Prom  (3.10)  and  the  first 
equation  in  (4.2),  assuming  for  simplicity  that  cw  =  0,  we  obtain  the 
rate  of  sublimation  in  the  final  formj 


L> « 


_  ,rl |  /  "i  Cl  .1.  _  'Vj  0 

•'«  ■  V  t  ~  1  <r„)  4- 


st:J  • 


where  cp^  is  the  specific  heat  for  a  constant  pressure  of  the  incident 


flow. 


The  mass  velocity  of  sublimation  will  be 


1/  r 

,  cT>j  Too) 

"if 

V  2.1/rr 

-<?(!)]  "iJ 

(4.4) 


The  temperature  profile  in  the  body  is  found  from  Formula  (3.13). 


5.  Numerical  Solution  of  the  Problem  for  a  -  0.7;  1  and  Sc  =*  1 

Solution  of  Problem  (3.4) -(3. 10)  in  the  general  form  is  associated 
with  numerous  calculations  owing  to  numerous  defining  parameters; 
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V  *v, 

M.  or  Sc,  — -  ,  .7-  (»«.  1 . A'-l).  Cm  (t-1 . A”).  *1. 

Vv  -v 


o.  e  T 

I  •  ” 

Pie" 


y»(0) 


,-/^c 

K 


(5.1) 


and  dimensionless  functions  PR  =  pr(S»  *).  Q  (0w.  ^(^VD.  *  (*«)/*«• 

We  therefore  made  certain  simplifications  which  reduce  the  number  of 
these  parameters. 

First  we  will  assume  that 

-l 


(5.2) 


This  corresponds  to  the  case  where  Le  Z  1  or  where  the  components  of 
the  mixture  have  similar  specific  heats  at  a  constant  pressure  (c^  3 
*  c  )  or  when  sublimation  is  weak  (c.  «  l)  .  Since  for  a  large 
number  of  gaseous  mixtures  these  assumptions  are  always  fulfilled  to 
a  certain  degree  simultaneously,  then  for  these  mixtures  assumption 
(5.2)  will  be  valid. 

Using  Formula  (2.3)  for  the  viscosity  of  the  mixture  and 


Sutherland* s  formula  for  the  viscosity  of  the  J-th  component,  we  obtain 

(f«.) M *  v  ev )  ’n> + (c?.lep>  1*0 + u  —  ’ 


—  !lP 
up  =  - -  1 

Hopo 


*} O'- 


v* 


>00 


(5.5) 


3  rv « 


where  SJ  is  the  Sutherland  constant  and  is  equal  to  107  for  T  -  97°K, 
180  for  T  -  970°K  and  825  for  T  =  4850°K  [15).  During  vaporization 
of  most  condensed  media  in  air  we  have  c  /c  >  1,  Vl/Ko  <1.  The 
ratio  qp  (cj)/q>  (c^o)  18  close  to  unity.  Therefore,  the  product  of 
the  first  three  factors  in  (5*5)  is  close  to  unity  and  for  certain 
mixtures  is  a  slightly  varying  function  of  the  composition  (this 
function  may  be  calculated  on  the  basis  of  the  formulas  of  (2). 
Therefore,  we  set  this  product  equal  to  some  average  constant  value 
and  designate  it  by 
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(5.*) 


A'** 


syto)  m 

(«fc) 


where  the  upper  bar  indicates  the  average  value  of  the  parameter 
between  its  value  at  the  outer  edge  of  the  boundary  layer  and  at  the 
sublimation  surface.  Obviously  the  parameter  K  is  unknown  until 
solution  of  the  problem.  Taking  into  account  also  that  c^  a  cp  we 
obtain 


l»?  «  A'-'O  (2),  0(=) 


+  (1  —  *»)  * 


(5.5) 


and  finally 


Pal?''1  =  So  +  (1  —  S„)  2. 


(5.6) 


We  note  that  this  approximation  will  effect  mainly  the  value  of 
<p"(0)  *  a.  But  since  g«(o)  and  z«(o)  ~  a1/3  it  follows  from  an 
asymptotic  solution  (6)  that  the  simplification  (5.6)  has  little  effect 
on  the  product  of  g*(0)  and  z » (  0)  ,  through  which  the  rate  of  sublima¬ 
tion  and  tenqperature  at  the  sublimation  front  are  determined  from  (5*7)  • 
After  these  assumptions.  System  (5.4)  is  rewritten  in  the  form 


+  n( ftp'  =  <j>'*  —  z0  —  (1  —  zj  z,  (5.7) 

(Oo-’zy  -f  iMpz'  =  0,  (O  Sc -'S'Y  +  rufg'  -  0.  (5.8) 

A  formal  solution  of  the  last  two  equations  in  terms  of  the  constants 
a  and  Sc  with  regard  to  (5.6)  and  (5.7)  will  be 


*  to) 


(n*  s.  »o.  a)  #•  «  ci>  (t):  Sc.  a) 

° 'Ty  “  «(oo;Sc,*^  a)’ 


(5.9) 


where 


r  * 

“  to;  t,  So.  a)  -  J<D-»  cxp(-  ntj?  to)  <J>“<to)<to.  8  «  9  (©>. 
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Determining  z « (0)  and  g* (0)  from  (5.9)  we  obtain  from  (5.7)  N  tran¬ 
scendental  equations 


—tiT'anaa  (oo;  a,  za,  a) 


(5.10) 


J‘oo  —  hn  o  <0  foo:  3,  za.  a)  **o 

(Tt)  =  "S  «  (oo;  Sc.  x*.  a)  e<*  —  eio ' 


(5.11) 


If  sublimation  proceeds  according  to  diffusion  kinetics  (f  >  0.1), 
then,  adding  the  first  Eq.  (3.88)  to  (5.10)  and  (5.H)  We  obtain  a 
system  of  N  +  1  equations  for  the  determining  of  N  +  1  unknowns: 
ct ,  Zq ,  c^  (i  *  N  ■■  l)« 

If  sublimation  is  nonequilibrium,  then,  adding  Eq.  (3. 10)  to  the 
system  (5.10),  (5.H)  for  N  =  2  we  obtain  a  system  of  three  equations 
for  the  determination  of  three  unknowns:  a,  z0  and  c0.  For  an  actual 
solution  of  systems  (5. 10),  (5.11)  and  (3.8)  it  is  necessary  to  compute 
the  function  o>(*;  r;  z0,  a)  .  For  Sc  =  1  and  n  =  1  the  system  (5.8) 
with  the  boundary  conditions 


<F(0)=a  <  0,  <p'(0)=;(0)=g(0)=  o.  <p'(oc)==;(oc)=g(oc)==l  .(,5.12) 


coincides  with  a  boundary-value  problem  flow  past  an  infinite  cylinder 
in  the  vicinity  of  the  stagnation  line  in  the  presence  of  slip  and 
injection  and  was  solved  numerically  by  Beckwith  [16]  for  a  *  0, 

-0.5,  -1,  for  z0  =  0,  0.5,  1  for  two  values  of  Prandtl  numbers  a  *  0.7 
and  1.  Numerical  integration  was  performed  for  4>(z)  «  1,  SjQ  **  0.2 
and  0.02.  In  addition,  for  o  -  0  the  results  of  numerical  integration 
for  the  axi symmetrical  case  (n  *  2)  with  a  *  Q.7  and  1  are  presented 
in  this  work  for  several  values  of  the  parameters  z0. 

The  results  of  Bechwith’s  calculations  [16]  for  the  case  of  a 
plane,  supplemented  by  numerical  integration  of  the  system  (5.8)  for 
the  axi symmetrical  case  where  a  -  -0.5,  -1,  o  -0.7,  1,  and  Sc  -  1  are 


presented  In  Tables  1-4.  Thus,  for  Sc  ■  1  and  a  «  0.7,  1  System  of 
Eqs.  (5.10),  (5.11) ,  (3.8)  or  (5.10),  (5.11),  (5.10)  together  with 
Tables  1-4  solve  the  given  problem  in  the  plane  and  axisymmetrical  case 

TABLE  1 

n  =  1  (plane  case) 


»(0)  -« 

1 

*•  • 

|  <<«)-• 

•  •1. 
O-l 

|  «-0.7 

O  - 1 

(j-t.tS 

0 

O.C'.SO 

0.C071 

•  - 

0 

0.05 

_ 

_ 

0.C89I 

0.4576 

0.5 

0.9548 

0.9562 

0.9109 

0.6894 

1 

1.2526 

1.2326 

1.2326 

1.2326 

0 

0.2007 

0.2966 

_ 

-0.5 

0.05 

_ 

'  _ 

0.3658 

0.1949 

0.5 

0.66C9 

0.6610 

0.6430 

0.6276 

1 

0.9602 

0.9692 

0.9692 

0.9692 

0 

0.0571 

0.0705 

—  ■ 

_ 

-1  i 

0.05 

— . 

_ 

0.1572 

— 

0.5 

0.4519 

0.4551 

0.445S 

0.4380 

’  1 

1 

0.7566 

0.7566 

0.7566 

0.7566 

TABLE  2  * 

n  *  1  (plane  case) 


O  (0)  ■ 

•• 

X’  (0),  ©  mm  I 

«*  (»)  I 

«o>  '  1 
® -!  | 

X*  (0) 

s'W 

x*(0) 

f«» 
t  (0) 

ir<o) 

IE* 

•  =  0.7 

o  -  0.7  { 

•  -  0.7. 

•  -  0.7, 

tj-mOIC 

•  -1 

(C.5) 

0 

0.5067 

0.43G2 

1.140 

1.127 

0 

0.03 

_ • 

— 

_ 

0.4661 

1.147 

0.2969 

1.160 

_ 

0.5 

0.5421 

0.4696 

1.146  ! 

0.4484 

1.150 

0.4317 

1.153 

1.139 

1 

0.5705 

0.495$ 

1.151  ; 

0.4958 

1.151 

0.4958 

1.151 

1.143 

0 

0.2031 

0.2103 

0.9439  | 

_ 

_ 

_ 

0.9818 

—".5 

0.05 

— 

_ 

—  ' 

0.2452 

0.9751 

0.1025 

0.8663 

_ 

0.5 

0.2580 

0.2600 

0.9S4G  ! 

0.2430 

0.9SI9 

0.2295 

0.9791 

1.015 

•1  •  • 

0.2950 

0.2934 

1.005  ! 

0.2934 

1.005 

0.2934 

1.005 

1.028 

0 

0.0211 

0.0496 

0.5565 

_ _ 

_ 

_ 

_ 

0.841  * 

—  1 

0.05 

_ 

_ 

_ 

0.0819 

0.C752 

0.0432 

_  • 

_ 

0.5 

0.0823 

fi.1112 

0.7419  1 

0.1010 

0.7326 

0.09:52 

0.7221 

0.892 

1 

0.1168 

0.1457 

0.8016  1 

0.1457 

0.SO16 

0.1457 

0.8O16 

0.918 
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TABLE  5 

n  *  2  (axisymmetrical  case) 


If  10)  —  a 

u 

»'  (0)  —  o, 

O-l 

*• 

.  •  »  0.7 

•  - 1 

•  —  0.7 

•j  —  0  02 

0 

0.8217 

0.7820 

0.05 

0.8C73 

0.5 

1 .070 

•  1 .055 

0.06 

— 

0.5820 

0 

0.0 

_ 

1.M0 

0:2 

— 

0.7483 

0.8 

— 

i.2«e 

0.5 

1.024 

0.9956 

* 

1.312 

1.312 

1 

1.312 

1.312 

0 

0.21S 

0.240 

— 

_ 

— 

A  C 

0.5 

0.551 

0.446 

— 

— 

— 

— U.*> 

0.8 

0.701 

0.702 

— 

— 

— 

1 

0.805 

O.S05 

— 

— i 

— 

—i 

0 

0.5 

0.281 

0.285 

I 

z 

1 

0.191 

0.491 

— 

— 

—  \ 

TABLE  ^ 


n  =  2  (axi symmetrical  case) 


6.  Solution  of  the  Problem  for  Any  Prandtl  and  Schmidt  Numbers 

In  order  to  obtain  a  solution  to  the  problem  for  any  values  of 
cr  and  Sc  we  will  evaluate  the  integrals  which  enter  into  system  (5.10)  4 
(5.11)  by  the  asymptotic  method.  Assuming  that  the  parameter  t  *  o 
(or  Sc)  is  sufficiently  large  and  that  parameter  a  is  sufficiently 
small  (a  -  1/r)  we  will  have 
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(6.1) 


where 


w  T,  co 

w(x;r.re.a)  =\  cxp(  wit  \<p  (r\)  <lr\)di\  exp(witai)  —  nil(r\))dr\. 


'=  Ol)  =  ^<T  (»l)  —  arj  =.  T]»  ^  <* ij\ 


£«  —  nin 


Co  — 


n*0 

»*  (s.4-  nan)  —  (1  —  *,)  *1 


“  ~  0  •  1  “  4!  ’  --  ~ - 5f 

_  (2  —  n)  n-  —  nV  (ru  -i-  nan)  -f  na  (1  -j-  c)  (1  —  s0)  x' 
c, - ;  j-i - — — - , 

«a(5n  —  8)  a:-rn[n‘a4-^if(4n— GJI+s^a*— nsa*(l+5+os)(l—  x,).x^ 

C*  71  — - —  , 


(6.2) 


«.  =  9*  (0).  v  -  2'(0). 

Changing  under  the  integral  in  (6.1)  to  the  variable  of  integra¬ 
tion  z  we  obtain  [17] 


®  m«0 


<»+«) 


(6.3) 


where 


■'--ra §«p (—1) r* 


V.  (m+1) 


From  Which,  using  (6.2),  we  obtain 


>7  *  T  Zn-'-nyr. 

(!a  =  1,  C.'i  =— c- - RTO, 


7  1  ntx'y~  nra  (s. -- wan)  ,  f:j  —  naa)5  na  (:u  --  naa)  —  (1  —  s.,)  xn 


•in  ■  10a*  20a  ’ 

*7  •  ft3!3!3  ,  n:TIa:fr5-'- nan)  ~nta  fa,  -j-  naa)*  ,  i»Ta|na(rg-;-na a) — (i — x.)x0]  , 

3  tj  *  <in  72a5  1  J5a  "7" 

,  So  (r3 -4- nan)5  7  fro •’an;  ,  ,  ,  .  ■ 

~  -TJTsT^ - "T^a—  l,,:t  (=»  +  ,laa)  “  V  ~  Z«)  z«]  ~ 

(2  —  n)  a5  —  n5a;  (r0 nan)  J- na  (I  J- s)  (i  —  s0) 


00a 


T  _  (.‘T'a*  on^rVfo -}- nanl  ,  5n*r5a3 (x#  -f-  nan)5 

4  2'i  725  ‘  725* 

«Vi! ! na  {s0  +  naa)  —  , '  — 1»)  z„ j 


24a 

_  M«Ta{s.-;.„an)»  ,  nta  (io  ~  naa)  [na  (:«  -f-  naa)  —  (1  —  s.)  *),]  t 

lO-bla3  »  lba*  "•* 

'•T3l<--,|),»i-,‘iai{-a-r»aa)-i-na(l-i-s)(1— -o).-,)  385  (sq  4-  naa)* 

*  72a  + 128-243  7* 
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li  (•*« /„  . _ \  ,4  *„  x_*, ,  i**  (**+»*«)— 0—4) 4i* 

; - itf-W"  l"*  (*•  T  na<*)  —  (1  —  2®)  "#)t - J5J75 - 

(2  —  n)  a*  —  n*a*  {/*  —  now)  +  ns  (1  -f  s)  (1  —  *,)  s' 


108  «* 


(*•  +  na  a)  — 


*a  (5»  —  8)  a*  +  a  I»*a*  +  (4*  —  6))  +  vrtt*  —  n*a* (1  +  o  +  s’)  0-4) rl 


(6.4) 


504  a 


Since  Zo  enters  Into  the  coefficients  32,  3s,  ...  then  expansion  of 
(6.3),  generally  speaking,  yields  a  transcendental  equation  for  the 
determination  of  o  in  terms  of  r,  z0  and  a.  But  since  the  main  con¬ 
tribution  to  the  value  of  the  function  co(»;  r,  Zo,  a)  is  introduced 
by  the  first  terms  of  the  series,  then  substituting  d2,  ...,  z&  ■ 

-  3  (ana/6)  1/ST“1  (j) ,  in  the  coefficinets,  sufficient  accuracy  can  be 
obtained.  Due  to  the  weak  dependence  of  the  magnltuded  of  a  on  a 
(see  Tables  1  and  3)  in  (6.3)  the  value  of  a  may  be  substituted  when 
a  -  0.7.  Calculations  according  to  Formula  (6.3)  indicate  that  for 
a  *  0  and  T  >  0.5  the  first  three  terms  yield  a  value  of  o>(»;  T,  z0,  0) 
with  an  error  less  than  Yf>  (see  Ref.  18)  .  When  a  «  -0.5  it  is  neces¬ 
sary  to  calculate  five  terms  of  the  series  (6.3)  in  order  to  obtain  an 
error  not  exceeding  ljt. 

By  expanding  (6.3)  and  limiting  ourselves  to  the  first  term,  we 
obtain 


<f), 


4 (°)  a («*:  g.  4. «)  .  fSc  \  I.  ,  V  »  // «  \  _ l_\  , 

:'(0)  u  (oo,  Sc,  4,  a)  \e/  |  *  L<*«  'j“*  Sc’-'*/"1" 

+  t  (jh  -  dv.)- *"  ^ -• Sc/’>] + "  *}• ' 


(6.5) 


Table  2  gives  for  comparison  the  values  of  g*(0)/z*(0)  as  cal¬ 
culated  according  to  Formula  (6.5).  When  |a|  £  0.5,  Formula  (6.5) 
gives  an  error  not  exceeding  4.5#  for  all  Zo.  When  a  ■  -1  the  error 
reaches  a  significant  value  and  in  Formula  (6.5)  it  is  necessary  to 
take  into  account  subsequent  terms.  From  the  system  (5.10),  (5.H) 
it  follows  that  the  case  of  a  ■  -1  corresponds  to  "low-enthalpy" 
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materials,  l.e.,  materials  having  (hj  -  ho)/j5»  1  and  a  high  vaporiza- 
tlon  rate.  Por  example,  when  a  «  1,  n  «  1  (h^  -  h0)/^«  <o(»,  e,  Zo*  -1) 
when  z0  *  0.5  and  ^7.7  when  z0  -  0.  Thus,  this  case  Is  of  little 
practical  interest  and  we  will  correct  formula  (6.5)  for  a  <  -0.5. 

Using  asymptotic  representation  of  (6.5),  Eq.  (5.11)  may  be 
written  in  the  form 


c-*;a 

TOT 


'(£■) 


/(«.*..  6.  SC) 


(6.6) 


where 


/  («.  V  s.)  -  1  +  0.500(if  [£  &-&}*?(£-£)- 
—  an  (cv*  —  Scv,)j- 


The  System  of  Eqs.  (5.10),  (6.6)  and  (5.8)  or  (5.10),  (6.6)  and  (3.10) 
determines  a  solution  to  the  problem  for  any  a  and  Sc.’ 


7.  Necessary  and  Sufficient  Conditions  for  Bolling  at  the 

Sublimation  Surface 

If  vaporization  is  equilibrium,  then,  as  was  pointed  out  In  4, 
boiling  on  the  sublimation  front  is  attained  when  there  is  Infinitely 
of  large  thermal  fluxes  from  the  gas. 

We  will  now  examine  the  conditions  of  reaching  the  boiling  point 
during  nonequilibrium  sublimation  (condition  3.10).  During  boiling, 
the  equilibrium  vapor  pressure  p(o)  becomes  equal  to  the  external 
pressure  p0o  and  the  temperature  at  the  sublimation  front  reaches  a 
maximum  value  T0  ■  T0o  for  a  given  external  pressure  p0o  at  the  stagna¬ 
tion  point.  Assuming  then  that  p0  *  Poo  and  T0  -  Too  in  (3.10)  and 
eliminating  the  parameters  a  and  Co  from  (5.10),  (5.H)  and  (3.10), 
we  obtain  the  necessary  and  sufficient  conditions  of  boiling  or  attain¬ 
ment  of  maximum  temperature  at  the  sublimation  front  for  a  binary 
mixture  (cw  -  0,  K  -  1,  Le  -  1)  In  the  form 
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i  7,«o “ ^ joo *'7  J \ 

!w1'tws;),s'u  ~-s#> — r 

nel  i 


(7.1) 


the  value  of  c0  which  enters  Into  the  parameter 

r.  «„  T. 

"0  P,  * 


K  ev  Tt  '*'»  +  %  t*~ 


ft*  ««_  ‘ 

00  Woo  40 


Woo 


.*0» 


being  determined  from  the  ratio 


e,  _  eP}lT<»~~T*)  (7.2) 

»-*“  WW3  * 

from  which  It  follows  that  there  are  always  finite  values  of  determin¬ 
ing  parameters,  satisfying  condition  (7.1)  at  which  boiling  ensures. 

In  fact,  when  the  parameter  (h^  -  hj0)/sf(Too)  la  increased  from  0  to 
*,  the  left-hand  part  of  Equality  (7*1)  increases  monotonically  from 
0  to  «  while  at  the  same  time  the  right-hand  part  decreases  monoton¬ 
ically  from  ocu(»;  c,  z0,  -(ns)”1)  >  0  to  ou>(w;  o>  z0,  0)  <  ».  Con¬ 
sequently  for  given  values  of  the  parameters  e,  M./M. ,  c  /c  and 

J  i  Pi  Pj 

# 

a  there  is  an  unique  value  of  the  parameter  (h^  -  hj0o)/^at  which 
boiling  seta  in. 

If  Le  /  1  then  the  necessary  and  sufficient  condition  of  boiling 
is  obtained  by  eliminating  parameter  a  from  the  system  of  two  trans¬ 
cendental  equations  (cw  «  0,  K  *  1) : 


V  *,</> f,  ,  ,lx~him/Se\  w~l  ,  1'  „  .1  . 

x  (TW)  L1  ~  x  (/’fc)  It;  1  '°’  zo'  Sc)J  ** au>  (°°;  zo>  °) 

t  r 

=..«.&»].+ 1-0, 


V* 


am 


(7.5) 


which  has  a  unique  solution. 

When  the  condition  (7.1)  Is  satisfied,  the  velocity  and  mass 
velocity  of  sublimation  are  found  from  (4.5)  and  (4.4)  respectively. 
When  Le  /  1  the  velocity  and  mass  velocity  of  sublimation  during 
boiling  may  be  obtained  according  to  the  formulas  (if  for  simplicity 
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we  set  I  si) 


ePj  —  ^#»)  J/j-  */•  "J 

l  {j'j  +  p,pr.‘e,  r-co  IQ  (Ol.)  -  Q  (0)  Wi  J  ' 

_ S(r«-ro»> _ *[+  Lc%  r* 

*  (^ce)  T  */-»  IQ  (Ou)  —  Q  (•))  A,i  p  J  * 
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